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Abstract The Large High Altitude Air Shower Observatory project is proposed to study high energy gamma 
ray astronomy ( 40 GeV-1 PeV ) and cosmic ray physics ( 20 TeV-1 EeV ). The wide field of view Cherenkov 
telescope array, as a component of the LHAASO project, will be used to study energy spectrum and compo- 
sitions of cosmic ray by measuring the total Cherenkov light generated by air showers and shower maximum 
depth. Two prototype telescopes have been in operation since 2008. The pointing accuracy of each telescope 
is crucial to the direction reconstruction of the primary particles. On the other hand the primary energy 
reconstruction relies on the shape of the Cherenkov image on the camera and the unrecorded photons due to 
the imperfect connections between photomultiplier tubes. UV bright stars are used as point-like objects to 
calibrate the pointing and to study the optical properties of the camera, the spot size and the fractions of 
unrecorded photons in the insensitive areas of the camera. 
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1 Introduction 

The Large High Altitude Air Shower Observatory 
( LHAASO ) [l| project aims to study 40 GeV-1 PeV 
gamma ray astronomy and 20 TeV-1 EeV cosmic ray 
physics at Yangbajing ( 4300 m a.s.l. ), Tibet, China, 
near the AS 7 and ARGO-YBJ experiments. The 
Wide Field of View ( FOV ) Cherenkov telescope Ar- 
ray ( WFCTA ), as a component of the LHAASO 
project, is designed to study cosmic ray energy spec- 
trum specie by specie by measuring the energy and 
X„, depth of each air shower. 

Two WFCTA prototype telescopes have been con- 
structed and placed nearby the ARGO-YBJ experi- 
mental hall. The two telescopes can be operated in 
both monocular and stereo modes, while coincident 
observation with the ARGO-YBJ detector is achieved 
off-line. 

Each telescope is made up of two main parts, the 



reflector and the camera. The reflector consists of 
20 spherical mirrors with a radius curvature R of 
4740±20 mm, corresponding to a total area of 4.7 
m 2 . The reflecting efficiency of the mirrors is about 
82% for light with wavelength larger than 300 nm. A 
camera is placed at the focal plane which is 2305 mm 
away from the reflector center to optimize the spot 
shape of a point-like object. The camera is composed 
of 256 flat hexagonal photomultipliers tubes ( PMTs ) 
each of which has a diameter of 40 mm, corresponding 
to a FOV of about 1° X 1°. The PMTs are arranged 
in 16 columns and 16 rows forming a total FOV of 
14° x 16° [2]. The maximum quantum efficiency of 
PMTs can reach 30% at 420 nm [H]. The signals of 
the PMTs are digitized by 50 MHz Flash Analog to 
Digital Converters ( FADCs ). The whole system is 
hosted in a shipping container with a dimension of 
2.5 mx2.3 mx3 m. The container is mounted on a 
standard dump-truck frame with a hydraulic lift that 



Received XX March XXXX 

* Supported by 100 Talents Programme of The Chinese Academy of Sciences, Knowledge Innovation Program of The Chinese 
Academy of Sciences ( H85451D0U2 ), National Natural Science Foundation of China ( 10975145 ). 
1) E-mail: llma@ihep.ac.cn 



No. X 



Lingling Ma et al: Geometry and optics calibration for WFCTA telescopes using star light 



2 



allows the container to be tilted in any elevation an- 
gles from to 60 degrees. The pointing direction of 
the telescope can be easily changed Q. 

The pointing accuracy and geometry properties 
of the telescopes are crucial to the arrival direction 
and X max reconstructions. The energy reconstruc- 
tion relies on the recorded Cherenkov image, how- 
ever, due to the imperfect physical junction between 
PMTs, some Cherenkov photons are unrecorded in 
the joints. In order to improve the accuracy of energy 
reconstruction the optical properties are studied. 

In this paper, we describe a method to calibrate 
the geometry and optical properties by using UV 
bright stars. While monitoring the UV Cherenkov 
light from air showers, the telescopes are also sensi- 
tive to the UV light from stars crossing the FOV of 
the telescopes. With their well known positions, or- 
ders of magnitude more accurate than the required 
resolution of WFCTA, and their point-like shape, the 
stars are ideal tools to test the pointing direction of 
each telescope. Using stars the optical properties of 
the telescope are studied. 

2 Night sky background and star sig- 
nal 

In addition to record the Cherenkov light from air 
showers, the camera also records the night sky back- 
ground light (NSB). When a bright star enters the 
FOV of the telescope, light from the star is added to 
the diffused NSB. The changes of the recorded NSB 
reflect the light sources crossing the FOV of the PMTs 
and the weather conditions. A star stays about 4 
minutes in the FOV of a PMT, during this time light 
from the star is added to the diffused NSB. While 
NSB change due to weather change behaves differ- 
ently, which usually lasts much longer than 4 min- 
utes and almost all PMTs are affected at the same 
time. This enables star signals to be discriminated 
from weather change. 

For each air shower event, the signal of Cherenkov 
light only lasts a few nanoseconds, while the trigger 
window lasts 18/iS, so telescopes record NSB in most 
of the trigger window. In order to reduce the fluctua- 
tion the recorded NSB is averaged in every 10s. Fig.Q] 
shows a typical NSB curve in one night recorded by 
a PMT. After subtracting the diffused background, 
peaks due to stars are clearly seen. The peak ampli- 
tude of a star light curve in a PMT depends on its 
UV brightness and its projected position to the PMT 
on the camera. In a typical night, many stars can be 
seen by each PMT. 



54739.6 54739.65 54739.7 54739.75 54739.8 54739.85 54739.9 



54739.6 54739.65 54739.7 54739.75 54739.8 54739.85 54739-jl 

Fig. 1. A typical NSB curve recorded by a PMT 
in one night before (upper) and after (bottom) 
subtracting the diffused NSB. 



3 Pointing direction 

Stars with well known positions and brightness are 
used as guiders of the telescopes. In our analysis, the 
TD 1 catalog is used, which has four different wave- 
length bands, 1565, 1965, 2365 and 2740 A respec- 
tively Since the WFCTA telescopes are sensitive 
in the near UV band, stars with flux in the 2740A 
band above 1 x 10 _11 erg/cm 2 / s are used. 

The pointing direction of each telescope can be 
changed through the container which encloses the 
whole telescope system. Only very rough pointing 
direction (about 200° in azimuth and 60° in eleva- 
tion) is known through it. Using the rough pointing 
direction, the time when the brightest star appears in 
the FOV of the telescope can be found through the 
brightest PMT. However, due to the large size of the 
PMT, the position of the PMT can not show accu- 
rately the position of the brightest star on the camera. 
So the weighted center position (x , y ) of the PMT 
and its neighbors within 2° is considered as the posi- 
tion of the star on the camera. When the star is in the 
middle of the camera in the horizontal direction, the 
star has the same azimuth angle with the telescope, 
while the elevation angle of the star is equal to the 
elevation of telescope plus the distance between the 
position of the star and the center of the camera. The 
obtained pointing direction of the telescope is more 
accurate than the one from the container. 

After getting the pointing direction of the tele- 
scope through the brightest star, orphan stars which 
have no surrounding stars within 2° are used to cor- 
rect the pointing direction of the telescope by using 
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iteration method. 

The accuracy of the pointing direction can be de- 
scribed by differences between positions of stars in 
the local coordinates and the obtained ones. Fig. [5] 
shows the distribution of the differences for one tele- 
scope. An accuracy better than 0.05° is obtained in 
less than 20 minutes with five stars in the FOV. 
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Fig. 2. The distributions of the differences be- 
tween positions of stars and the obtained ones 
in azimuth (upper) and elevation (bottom) . A 
Gaussian fit is superposed to each distribu- 
tion. 



Fig. |3] shows the pointing of one telescope chang- 
ing with time in two observation periods, i.e. two 
months. Both azimuth and elevation are very stable. 
The big change of the elevation was due to lowering 
and elevating the container between the two periods. 
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Fig. 3. The pointing directions of one telescope 
versus time in two months (Up: azimuth, Bot- 
tom: Elevation). Each point represents one 
day. 



4 Camera geometry 

The camera geometry calibration is done after 
pointing correction of each telescope. The calibra- 
tion includes the following four parameters. The first 
one ( Pi ) is a scaling of the tubes away from the cen- 
ter of the tube cluster. To first order P x corrects for 
deviations in the radius curvature of a mirror and for 
changes in the effective camera-mirror distance due to 
the treatment of the flat camera as a curved surface. 
The second parameter ( P 2 ) describes the rotation 
angle of the camera around the mirror axis. The last 
two parameters ( P A , P 4 ) indicate the offsets of the 
shift in the position of the entire camera with respect 
to the mirror axis. After the four parameters correc- 
tions the positions (x c ,t> y c ,t) of stars on the camera 
at time t are modified by Eq. (fTJ) and Eq. ([2]). 

x' ct = (l+pi)(x Cit cos(p 2 )-y Cit sin(p 2 ))+p3 (1) 

y'ct = (1 +Pi)(xc,tsin(p 2 ) + y c , t cos(p 2 )) +p4 (2) 

The x' c t and y' ct in Eq. (fTJ and Eq. © are the 
stars' coordinates after geometry correction. The 
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four parameters can be obtained by the least squared 
method. The x 2 is shown in Eq. <J3]) . 



X 



EE 



(x' t X startt ) 2 + (y' ct Ustar,tY 



(3) 



In Eq. ([3]), the x ata r,t, Dstar,t are the stars' exact 
positions on the camera at time t. The a x and a y are 
the RMS of the x' ct and y' c t , respectively. By min- 
imizing the x 2 , the values of P ll P 2 , P3 and P4 are 
found to be -1.5%, 0.6°, -0.05° and -0.05°. These 
parameters are used in the detector simulation and 
data reconstruction. 

5 Spot size and variations of the num- 
ber of observed photons 

Photons from any given direction form a quasi- 
Gaussian-shaped spot on the camera, which is in- 
fected by imperfection of the mirrors surface. The 
spot size depends on the angular distance to the op- 
tical axis. The larger the angular distance to the 
optical axis, the larger the spot size and the more it 
deviates from a Gaussian shape due to coma of the 
image. The spot size as an important parameter that 
affects the Cherenkov images of air showers is taken 
into account in the ray tracing procedure in both data 
analysis and detector simulation. 





Fig. 4. One example of light curves in X (left) 
and Y (right) directions when a star passes 
through a PMT's FOV. 

Bright stars can be considered as perfect point 
sources. The light from a bright star is collected by 
mirrors and projected to the camera, forming a light 
spot. The camera records the light spot in poor res- 
olution due to the large pixel size. If a PMT is on 
the track of a bright star, the PMT gets brighter 
and brighter as the star goes nearer, dimmer and 
dimmer as the star leaves. If the spot size is much 
smaller than the track length of the star in the FOV 
of the PMT, the light curve recorded by the PMT will 



be rectangular-shaped, with the width approximately 
equal to the track length, while if the spot size is much 
larger, the light curve also will be rectangular-shaped, 
with the width equal to the spot size. In our case, the 
spot size is similar to the pixel size, the light curve 
becomes quasi-Gaussian-shaped. 

Fig. @] shows one example. The spot size can be 
estimated by fitting the light curve accordingly. In 
order to get rid of the affects of nearby stars, only or- 
phan stars, which have no nearby stars within 2°, are 
used. Fig. [5] shows the variation of the spot size with 
angular distance to mirror axis. The spots become 
larger from the center of the camera to the edge. The 
spot size obtained from bright stars is used to improve 
the energy and arrival direction reconstruction. 
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Fig. 5. Spot size varies with angular distance 
to mirror axis. 

Photons which fall in the insensitive areas of the 
camera are never recorded. In the data analysis 
and detector simulation, the unrecorded photons have 
been taken into account as part of the ray tracing 
procedure. This is important in shower energy esti- 
mation and tested using the stars crossing the field 
of view of the telescopes. When a bright star with a 
constant flux passes through the FOV of the camera, 
the number of recorded photons varies due to the dif- 
ferent positions of the star on the camera. When the 
star is near to the center of a PMT, most photons 
fall in the sensitive areas of the camera, while when 
the star is near the conjunction of two PMTs, most 
photons fall in the insensitive areas of the camera, so 
the variations of the recorded photons from the star 
on its track can be observed which is shown in the 
Fig.lH The variations are also effected by the weather 
conditions, so several clear nights data is used to get 
an average behavior of it. In the Fig. [6j the simu- 
lation of the variations along the track of the star is 
also shown which is consistent with the observed one. 
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This demonstrates that the ray tracing simulation 
correctly copes with the fraction of the unrecorded 
photons in the insensitive areas of the camera. 
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Fig. 6. The back squares show the variations 
of the observed photons on the track of a star, 
while the black curve shows the simulated one. 

6 Conclusions 

Bright stars are used as a guider to calibrate the 
pointing direction of each telescope and the geome- 
try of the camera. The pointing accuracy obtained 
through bright stars is better than 0.05°. The long 



term stability of the pointing direction of the tele- 
scope is also monitored by bright stars. Moreover, as 
point sources, the bright stars are also used to study 
the spot shape. The spot size becomes larger from 
the center to the edge of the camera. Besides the 
spot size, the fraction of the unrecorded photons in 
the insensitive areas between PMTs is compared with 
the simulated one, and they are consistent with each 
other. The errors caused by this effect are well under- 
stood and under control in the energy reconstruction. 

The pointing direction and the correction of the 
geometry and optics of the telescopes are used in the 
simulation and data analysis to improve the recon- 
structions of the energy and arrival direction of the 
air shower. 
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